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Abstract 

Spin-orbit coupling is often described in an approach known as "the MacDonald 
torque", which has long become the textbook standard due to its apparent simplicity. 
Within this method, a concise expression for the additional tidal potential, derived by 
MacDonald (1964; Rev. Geophys. 2, 467 - 541), is combined with a convenient as- 
sumption that the quality factor Q is frequency- independent (or, equivalently, that the 
geometric lag angle is constant in time). This makes the treatment unphysical because 
MacDonald's derivation of the said formula was, very implicitly, based on keeping the 
time lag frequency-independent, which is equivalent to setting Q to scale as the inverse 
tidal frequency. This contradiction requires the entire MacDonald treatment of both 
non-resonant and resonant rotation to be rewritten. 

The non-resonant case was reconsidered by Efroimsky & Williams (2009; Cel.Mech. & 
Dyn.Astr. 104, 257 - 289), in application to spin modes distant from the major com- 
mensurabilities. In the current paper, we continue this work by introducing the necessary 
alterations into the MacDonald-torque-based model of falling into a 1-to-l resonance. 
(The original version of this model was offered by Goldreich 1966; A J 71, 1 - 7.) 

Although the MacDonald torque, both in its original formulation and in its corrected 
version, is incompatible with realistic rheologies of minerals and mantles, it remains a 
useful toy model, which enables one to obtain, in some situations, qualitatively meaningful 
results without resorting to the more rigorous (and complicated) theory of Darwin and 
Kaula. 
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We first address this simplified model in application to an oblate primary body, with 
tides raised on it by an orbiting zero-inclination secondary. (Here the role of the tidally- 
perturbed primary can be played by a satellite, the perturbing secondary being its host 
planet. A planet may as well be the perturbed primary, its host star acting as the tide- 
raising secondary.) We then extend the model to a triaxial primary body experiencing 
both a tidal and a permanent-figure torque exerted by an orbiting secondary. We consider 
the effect of the triaxiality on both circulating and librating rotation near the synchronous 
state. Circulating rotation may evolve toward the libration region or toward a spin faster 
than synchronous (the so-called pseudosynchronous spin). Which behaviour depends on 
the orbit eccentricity, the triaxial figure of the primary, and the mass ratio of the secondary 
and primary bodies. The spin evolution will always stall for the oblate case. For libration 
with a small amplitude, expressions are derived for the libration frequency, damping rate, 
and average orientation. 

Importantly, the stability of pseudosynchronous spin hinges upon the dissipation model 
employed. Makarov and Efroimsky (2012; arXiv:1209.1616) have found that a more re- 
alistic tidal dissipation model than the corrected MacDonald torque makes pseudosyn- 
chronous spin unstable. Besides, for a sufficiently large triaxiality, pseudosynchronism is 
impossible, no matter what dissipation model is used. 

1 Motivation 

Bodily tides in a near-spherical homogeneous primary perturbed by a point-like secondary are 
described by the theory developed mainly by Darwin (1879, 1880) and Kaula (1966). Sometimes 
this theory is referred to as the Darwin torque. Based on a Fourier-like expansions of the 
perturbing potential and of the tidally-induced potential of the disturbed primary, their theory 
permits an arbitrary rheology of the primary. 

Sometimes a much simpler empirical model, offered by MacDonald (1964) and often named 
as the MacDonald torque, is used in the literature for obtaining very approximate but still 
qualitatively reasonable description of tidal evolution. This model can be derived from the 
Darwin-Kaula theory, under several simplifying assumptions. Among these is the assumption 
that the tidal quality factor Q of the perturbed primary should scale as the inverse of the 
tidal frequency. Being incompatible with the rheologies of actual minerals, this key assumption 
prohibits the use of the MacDonald model in long-term orbital calculations. 

Despite the unrealistic rheology instilled into the MacDonald model, the model remains 
a lab which permits one to gain qualitative understanding of tidal evolution (rotational and 
orbital), without resorting to the lengthy calculations required in the accurate Darwin-Kaula 
approach. It should be noted however that employment of the MacDonald model needs some 
care. Historically, the orbital calculations performed with aid of this model by one of its cre- 
ators, MacDonald (1964), contained an inherent contradiction. MacDonald began his paper 
with deriving a concise expression for the additional tidal potential, and then combined this ex- 
pression with a convenient assumption that the geometric lag angle is constant (or, equivalently, 
that the tidal quality factor is a frequency independent constant). This made the treatment 
unphysical because MacDonald's derivation of the said formula was, very implicitly, based on 
keeping the time lag frequency-independent, which is equivalent to setting Q to scale as the 
inverse tidal frequency. This contradiction made his theory inconsistent. The said oversight 
has then been repeated many a time in the literature (Kaula 1968, eqn. 4.5.37; Murray & 
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Dermott 1999, eqn. 5.14). In particular, the MacDonald method was used by Goldreich 
(1966) in his theory of dynamical evolution near the 1:1 spin-orbit resonance. We reconsider 
this theory, employing the corrected version of the MacDonald torque, i.e., setting the quality 
factor to scale as inverse frequency. 

2 Linear bodily tide in a near-spherical primary 

Consider a near-spherical primary of radius R and a secondary of mass M* ec located at r* = 
(r*, 0*, A*), where r* > R. The tidal potential created by the secondary alters the primary's 
shape and, as a result, its potential. For linear tides, the amendment to the primary's exterior 
potential is known (e.g., Efroimsky & Williams 2009, Efroimsky 2012a,b) to read as 

U(f) = -GMl^kt l+1 l+1 ^) y ^ T (2-^)^(sin0)P ;m (sm0*) cosm(A-A*) , (1) 
i =2 r r* m=0 \ L + m )- 

5ij being the Kronecker delta, G = 6.7 x 10~ n m 3 kg~ 1 s -2 being Newton's gravity constant, 
and 7 being the angle between the vectors r* and r pointing from the primary's centre. As 
agreed above, r * = (r* , <fi* , A*) denotes the position of the perturber, while r = (r , <fi , A) is 
an exterior point, at a radius r > R, where the tidal potential amendment U(r) is measured. 
The longitudes A, A* are reckoned from a fixed meridian on the primary, while the latitudes 
0, 0* are reckoned from its equator. Indices I and m are traditionally referred to as the degree 
and order, accordingly. The associated Legendre functions Pi m (x) (termed associated Legendre 
polynomials when their argument is sine or cosine of some angle) are introduced as in Kaula 
(1966, 1968), and may be called unnormalised associated Legendre functions, to distinguish 
them from their normalised counterparts!! The Love numbers ki can be calculated from the 
geophysical properties of the primary. 

A different formula for the tidal-response-generated change in the potential was suggested 
by Kaula (1961, 1964). Kaula devised a method of switching variables, from the spherical 
coordinates to the orbital elements (a*, e*, i*, Q*, cu*, Ai*) and (a, e, i, Q, to, M. ) of the 
secondaries located at r * and r. The goal was to explore how a tide-raising secondary at r * 
acts on a secondary at r through the medium of the tides it creates on their mutual primary. 

The development enabled Kaula to process ([1]) into a series, which was a disguised form of 
a Fourier expansion of the tide. Interestingly, Kaula himself never referred to that expansion 

1 Due to an error in our translation from German, we mis-assumed in our previous papers Efroimsky & 
Williams (2009) and Efroimsky (2012a) that Gerstenkorn (1955) had based his development on a constant- 
Q model. Therefore we stated that his theory contained the same genuine inconsistency as the theory by 
MacDonald (1964). Accurate translation of the work by Gerstenkorn (1955) has shown that his method was 
based on a constant-time-lag model. Therefore we retract our statement about Gerstenkorn's approach sharing 
the inconsistency of MacDonald's theory We also thank Hauke Hussmann and Peter Noerdlinger for their kind 
help in translating excerpts from Gerstenkorn's work. 

2 The Legendre polynomials may be defined through the Rodriguez formula Pi{x) = —f-r -j-j ( x 2 — l)' . 
In most literature, the associated Legendre functions are introduced, for a nonnegative m , as 

P lm (x) = (1 - x 2 ) m/2 ^mx) and PT(x) = (-1)'" (l - ^ , 

so that Pim(x) = (—l) rn P[ n (x). The above definition agrees with the one offered by Abramowitz & Stegun 
(1972, p. 332). A different convention is accepted in those books (e.g., Arfken & Weber 1995, p. 623) where 
P; m (x) lacks the (— l) m factor and thus coincides with Pi m {x) . 
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as a Fourier series, nor did he ever write down explicitly the expressions for the Fourier modes. 
At the same time, the way in which Kaula introduced the phase lags indicates that he was 
aware of how the modes were expressed via the perturber's orbital elements and the primary's 
rotation rate. The original works by Kaula (1961, 1964) were written in a terse manner, with 
many technicalities omitted. A comprehensive elucidation of his approach can be found in 
the Efroimsky & Makarov (2012). Referring the reader to that paper for details, here we cite 
only the resulting formula for the secular part of U , in the special case when the tide-raising 
secondary itself experiences perturbation from the tide it creates on the primary (so r* — r): 

U {sec \r) = 



21+1 , , 
E - E 77-S ( 2 - W E F U0 E k ^ Pq ) COSe ; (^ mp9 ) , (2) 



G M SP . r . / R\ Tr^il-m, 

ry (2 — d 0m ) 

p=0 



where I, m, p, q are integers, Fi mp (i) are the inclination functions (Gooding & Wagner 2008), 

Gipq(e) are the eccentricity polynomials coinciding with the Hansen coefficients X^_ 2p+ ^ (e) , 
while the superscript "sec" means: secular. 

The dynamical Love numbers ki and the phase lags e; are functions of the Fourier modes 

ui mpq = (I -2p) Cj + (I -2p + q) M + m (Cl - 6) , (3) 

9 being the primary's sidereal angle, and 9 being its angular velocity. While the tidal modes 
(J3J) can be of either sign, the physical forcing frequencies 

Xlmpq = I Ulmpq \ = | (I - 2p) U + (I - 2p + q) M + 171 (fi - 9) \ (4) 

at which the stress and strain oscillate in the primary are positive definite. 

A partial sum of series (j2J), with |/|, \q\, \j\ < 2 , was offered earlier by Darwin (1879). An 
explanation of Darwin's method in modern notation can be found in Ferraz-Mello, Rodriguez 
& Hussmann (2008)0 

The power of the Darwin-Kaula approach lies in its compatibility with any rheology, i.e., 
with an arbitrary form of the mode-dependence of the products ki cos e\ . It can be demon- 
strated that, for a homogeneous near-spherical primary, the functional form of the mode de- 
pendence of the product k\ cose/ mpq is defined by index / solely: ki(oji mpq ) cos ei(ui mpq ) , the 
other three indices being attributed to the tidal mode. One can assume that this product de- 
pends not on the tidal mode U[ mpq but on the positive definite physical frequency Ximpq ■ This 
however will require some care in derivation of the tidal torque from the above expressions for 
the potential - see Efroimsky (2012a,b) for details. 

We had to write down the secular part of the Kaula expansion of tide, because we shall use 
it as a benchmark wherewith to compare the empirical expression by MacDonald (1964), which 
we shall derive below in an accurate manner. 

Lacking the ability to accommodate an arbitrary rheology, the MacDonald approach (or 
the MacDonald torque) produces, after a necessary correction, a simple method which can 
sometimes be employed for getting qualitative understanding of the picture. 



In Ferraz-Mello, Rodriguez & Hussmann (2008), the meaning of notations r and r * is opposite to ours. 



4 



3 Tidal torque 



Consider a secondary body of mass M sec , located relative to its primary at r = (r, A, <fi) , 
where <fi is the latitude, and A denotes the longitude reckoned from a meridian fixed on the 
primary. Let U stand for the tidal-response amendment to the primary's potential. This 
amendment can be generated either by this secondary itself or by some other secondary of 
mass M* ec located at r* = (r*, A*, 0*) . In either case, the primary's tidal response to the 
gravity of the secondary will render a tidal force and torque acting on the secondary of mass 
M sec • The torque's component perpendicular to the equator of the primary will be given by: 

T x = - M m ?W. . (5) 

We would reiterate that (jSJ) is a component of the torque wherewith the primary acts on the 
secondary of mass M sec . 

Then its negative will be the appropriate (i.e., orthogonal to the primary's equator) com- 
ponent of the torque wherewith the secondary acts on the primary: 

T z {f) = — T z = M sec ^p- . (6) 

Derivation of formulae (jS - E} is presented in Appendix |A] These expressions are convenient 
when the tidal-response potential amendment U is expressed through the spherical coordinates 
r , A , and r * , A * , * , as in formula flTJ . 

Whenever the tidal response is expressed as a function of the orbital elements of the 
secondary and the sidereal angle 9 of the primary, it is practical to write down that the 
perpendicular-to-equator component of the torque acting on the primary as 

T,(f) = - M sec ^p- , (7) 

9 standing for the primary's sidereal angle (Efroimsky 2012a,b). 

We prefer to employ the terms primary and secondary rather than planet and satellite. 
This choice of terms is dictated by our intention to apply the below- developed machinery to 
research of tidal dissipation and spin evolution of a satellite. In this setting, the satellite 
is effectively playing the role of a tidally-distorted primary, its host planet acting as a tide- 
generating secondary. Whenever the method is applied to exploring the problem of planet 
despinning, the planet is understood as a primary, the host star being a tide-raising secondary. 



4 MacDonald (1964) 

When it comes to taking dissipation into account, expression ([I]) turns out to be far more 
restrictive than (j2J), in that ([T|) becomes applicable to a very specific rheological model. This 
happens because a straightforward^ option of instilling the delay into ([1]) is to replace in this 

4 A more accurate treatment, which cannot be employed within the MacDonald model but is implementable 
within the Darwin-Kaula approach, is to expand the tide-raising potential W and the tidal-response potential 
change U into Fourier modes uii mpq , and then to introduce the Love numbers ki mpq = ki(u>i mpq ) , phase lags 
timpq = £i(wi m pq) , and time lags Ati mpq = Ati(oji mpq ) . This formalism (explained in detail in Efroimsky 
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expression the perturber's coordinates r*(t) , <$>*(t) , X*(t) with their delayed values, r*(t — 
At) , <p*(t — At) , \*(t — At) . For example, instead of cosm(A — A*) we should employ 

/ (delayed) \ , . 

/ * \ / * . * \ 

cos m A — m A = cos I m A — ( mA - m A Af ) , (8) 



insertion whereof into ([I]) yields: 

°\ r ) , J+i , • J-M Z-e c; I mil ^ 



-5 Om )Pi m (sin0(t))P Zm (sin0*(t-At)) cosm(A - A* + A* At) . (9) 
Expressing the longitude via the true anomaly v , 

\ = -e + n + uj + u + 0(i 2 ) = -6 + Vt + u + M + 2e mxM + 0(e 2 ) + 0(i 2 ) 
and neglecting the apsidal and nodal precessions, we obtain: 

cos y ( m A - mA* ) + mA* At^j = cos ( ( m A - mA* ) + m (z>* - 0*) At + O (i* 2 ) 
or, in terms of the mean anomaly: 

cos ( (mA — mA ) + m \ * At 



cos ((mA-mA*) +m(n* - d*)At + 2me*n*At cosM*+ 0(e* 2 ) + O (i* 2 ) + 0(z 2 )) 
This enables us to write down the potential as 

no / 

f? ( 1 \ I 

17(f) = -GJiCc E*« a ^ Et7tSt ( 2 - 

1=2 r(t) r(t-At) m=G V^ m )- 

Pi m (sin c(>*(t- At)) cos(m(A-A*)+ m(v*-6*)At + 0(i 2 ) + 0(i* 2 )^) , (13) 

where P; m (sin 4>(t) )P; m (sin 4>*(t-At) ) may be replaced, in the order of 0(i 2 )+0(i* 2 )+0(ii*) , 
with Pim(0)Pim{0) ■ A further simplification can be achieved by taking into account the I = 2 
contribution only Omitting the A-independent term with m = (in the expression for the 



2012a, b) enables one to express the so-introduced Love numbers through the rheological properties of the 
primary's mantle, and thereby to model adequately the frequency-dependence of these Love numbers. 

An intermediate, purely empirical, option would be to introduce "Love numbers" fc; m , as if they were 
functions of both the degree I and order m . This idea is implemented in the IERS Conventions on the Earth 
rotation (Petit & Luzum 2011). In the LLR (Lunar Laser Ranging) integration software, tides in the Earth are 
parameterised by k\ m and Ati m with I = 2 and m = 0, 1, 2 (Standish and Williams 2012). 



6 



torque, m will become a multiplier after the differentiation of U(r) with respect to A), and 
omitting the m = 1 term (as 72i(0) = ), we arrive at the expression 

5 

U(r) = - - G ^e C k2R r x + Q( . 2) + Q( .* 2) + Q( C os f ( 2 A - 2A* ) 
4 r(t) r*(t - At) v 

+ 2 (z>* - 0*) At + 0(i 2 ) + 0(r 2 ) ) . (14) 

In the special case when the tide-raising satellite is the same body as the tidally-perturbed 
one (i.e., when r(t) = r*(t) , M sec = M* ec , and A = A* ), this expression happens to coincide, 
in the leading order of i , and i* , with the potential employed by MacDonald (1964). Thus 
we have reproduced his empirical approach, by starting with the rigorous formula ([!]), and by 
performing the following sequence of approximations: 

• First, when accommodating dissipation, we set all the time delays to be equal, for all the 
physical frequencies Xi mpq — I ^impq I involved in the tide: 

&t lmpq = At( Xlmpg ) = At . (15) 
This point is explained in great detail in Efroimsky & Makarov (2012). 

• Second, we assume the smallness of the inclinations and lag, through the neglect of the 
relative errors 0(i 2 ) , 0(i* 2 ) , and 0(ii*) . 



• Third, we truncate the series by leaving only the I — m — 2 term. 



These three steps take us from (JTJ) to the approximation 

17(f) « - - A ° f sec k2 R " , cos ( 2 A - 2A* + e ) , (16a) 
4 r(tfr*(t-Atf V ' 

With the tide-raising secondary set to coincide with the one perturbed by the tides on the 
primary (so r{t) = r*(t) , M sec = M* ec , and A = A* ), the above expression assumes the form 
of 

5 

U(r) w - - G M sec k 2 , R — cose . (16b) 

v ' 4 r(tf r(t -At) K ' 

In the denominator, r(t — At) can be replaced, in the order of 0(e*n* At) , with r*(t) . With 
this simplification implemented, we end up with 

5 

17(f) « - ? G M sec k 2 cose . (16c) 



5 Prom r = a(l — e 2 )/(l + e cosv) and dv/dM = (1 + e cos^) 2 /(l — e 2 ) 3 / 2 it is straightforward that 

a , \ , a \ a e (1 — e 2 ) . . ae sin v A _ / . ..on 

Ar = r(t) - r(t - At) = -- — i '— sin v Au + O (e 2 ) = - — - nAt + 0(e(n At) 2 ) . 

(1 + e cosv) z (1 — e z ) L ' z 

Thus our replacement of r*(< — At) with r*(t) = r(i) entails a relative error of order O(enAt) . In expressions 
(O, p^)) . ([13]), and P^)l . the absolute error will be of the same order, for A ^ A* . However for A = A* the 
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Expressions (|T4|) and ( TT6]) contain the longitudinal lag 



e = m \ At = 2(z>* - 9) At + 0{i 2 ) , (17a) 

v* and 9 being the true anomaly of the perturber and the sidereal angle of the primary. In 
the special case ( I16bl - [T6"cj) . when the tide-generating secondary and the secondary disturbed 
by the tide on the primary are one and the same body, the asterisks may be dropped: 

e = m A At = 2 (z> - 9) At + 0{i 2 ) . (17b) 

The spin rate 9 is a slow variable, in that it may be assumed constant over one orbiting 
cycle. The true anomaly is a fast variable. For a nonvanishing eccentricity, v too is a fast 
variable, and so is the lag e . This means that we should take into account these two quantities' 
variations over an orbital period. 

Expression f|16bj) coincides, up to an irrelevant constant J§ with the leading term of the 



appropriate formula from MacDonald (1964, eqn 21). To appreciate this fact, notice that, up 
to 0(i 2 ) , the absolute value of the longitudinal lag e is the double of the geometrical angle 
5 = \{9 — i>) At\ subtended at the primary's centre between the directions to the secondary and 
to the bulgeQ provided At is postulated to be the same for all tidal modes. 

The analogy between the MacDonald theory and that of Darwin and Kaula can be traced 
also by starting from the series ($Z§ . Consider the case of the tidally perturbed secondary 
coinciding with the tide-raising one, so A = A* , and all the orbital variables are identical to 
their counterparts with an asterisk. It will then be easy to notice that, formally (just formally), 
expression (Till) mimics the principal term of the series (j2]) , provided in this term the multiplier 



G200 * s replaced with unity, and the principal phase lag £ 22QQ = 2 (n — 9) At 22Qo is replaced 
with the longitudinal lag ( TTT|) . This way, within the MacDonald formalism, the longitudinal lag 
(flTl) is playing the role of an instantaneous phase lag associated with double the instantaneous 
synodic frequency 

X = 2 \u - 9\ , (18) 

which is, up to 0(i 2 ) , the double of the angular velocity wherewith the point located under 
the secondary (with the same latitude and longitude) is moving over the surface of the primary. 

To extend further the analogy between the MacDonald and Darwin-Kaula models, one can 
define an auxiliary quantity 

U Q" = —^TT ( 19 ) 
sin e 



absolute error will become 0(enQ~ l At) , since sine is of the same order as the inverse quality factor Q . In 
our estimates of errors, it is irrelevant whether we define Q as that appropriate to the principal tidal mode or 
via formula (|19j) below; so we simply use the generic notation Q. 

As wc explained in Efroimsky & Williams (2009), after averaging over one revolution of a nonrcsonant 
secondary about the primary, the absolute error reduces to O (e 2 n 2 Q _2 (Ai) 2 ) . In a resonant case, though, we 
cannot enjoy this reduction, because in this case the averaging procedure looks different as the torque changes 
its sign over a period of the moon's revolution. For this reason, the absolute error remains 0(enQ~ 1 At) . 

6 The additional tidal potential U , as given by equation (21) in the work by MacDonald (1964), fails to 
vanish in the limit of zero geometric lag S . This minor irregularity, though, does not influence MacDonald's 
calculation of the tidal torque. 

7 Be mindful that the double of the geometrical angle is not equal to the absolute value of e 2 200 = 2(n — 

Q) ^2200 • 
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and derive from ([IT]) and (fT9l) that, in the leading order of e and z , this quantity satisfies 



V = — TT • (20a) 
X At 

A popular fallacy would then be to interpret (I20ap as a rheological scaling law Q ~ with 
a = — 1 . That this interpretation is generally incorrect follows from the fact that the quantity 
U Q" , defined through ( !T9|) . is not obliged to coincide with the quality factor]^ 

To sidestep these difficulties, it would be safer to write the constant-time-lag rheological 
law, for small lags, simply as 

|e| = xAt . (20b) 

Historically, MacDonald (1964) arrived at his model via empirical reasoning. He certainly 
realised that the model was applicable to low inclinations only. At the same time, this author 
failed to notice that the model also implied the frequency-independence of the time lag. In fact, 
this frequency- independence, (j!5J) . is necessary to derive the MacDonald model ffTB"]) from the 
generic expression ([1]) for the tidal amendment to the primary's potential. This way, equality 
( TT5l) . is a priori instilled into the model. In other words, the MacDonald model of tides includes 
in itself the rheological scaling law ( I20bl) with a constant At . Unaware of this circumstance, 
MacDonald (1964) set the angular lag to be a frequency- independent constant, an assertion 
equivalent to the time lag scaling as inverse tidal frequency. However, as we just saw above, a 
consistent derivation of the MacDonald tidal model requires that the time lag be set frequency- 
independent |^] Thus, to be consistent, the MacDonald method must be corrected by applying 



8 Within the Darwin-Kaula theory, for each tidal mode uJi mpq , we introduce the phase lag as eimpq = 
^impq ^timpq ■ Then we introduce the appropriate quality factor Qimpq = Qiximpq) = Q{ \^impq 

| ) via the 

expression for the one-cycle energy loss: 

^^cycleyXlmpq ) — ^ 

Finally, using the fact that Xi mpq = l^impq I ^ s ^ ne frequency of a sinusoidal load, we prove that the afore 
introduced e; mp<J and Qi mpq are interconnected as 1/Q = sine, if E pea k denotes the maximal energy, or in a 
more complex way, if E pea k stands for the maximal work (Efroimsky 2012a, b). 

Within the MacDonald method, original or corrected as (fT5)) . it is not a priori clear if the overall peak 
work (or the overall peak energy stored) and the overall energy loss over a cycle are interconnected via the 
auxiliary quantity "Q" in exactly the same manner as E pea k(xi m pq) an( i ^-^cycieiXimpq) are interconnected 
by QiXimpq) m the above expression. (Recall that the total cycle of the tidal load is, generally, nonsinusoidal.) 
Whenever we can prove that 

r, j-p(overall) 
* ^(overall) _ ^ peak 

cycle ' 



our "Q" can be spared of the quotation marks and can be called instantaneous quality factor, while (|20a[) 
can be treated as a reasonable approximation to scaling law Q ~ \ a with a = —1 . However, this should be 
justified in each particular case, not taken for granted. 

9 In application to a non-resonant setting, a constant- At approach was taken yet by Darwin (1879). Later, 
MacDonald (1964) abandoned this method in favour of a constant-geometric lag calculation. Soon afterwards, 
though, Singer (1968) advocated for reinstallment of the constant-At method. Doing so, he was motivated by 
an apparent paradox in MacDonald's treatment. As explained by Efroimsky & Williams (2009), the paradox is 
nonexistent. Nonetheless, the work by Singer (1968) was fruitful at the time, as it renewed the interest in the 
constant-At treatment. The full might of the model was revealed by Mignard (1979, 1980, 1981) who used it 
to develop closed expressions for the tidal force and torque. 
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( JT5|) or, equivalently, ( 120b[) . In greater detail, the necessity of this amendment is considered in 
Efroimsky & Makarov (2012). 

Even after the model is combined with rheology (fT5|) . predictions of such a theory are of 
limited use. The problem is that this rheology is radically different from the actual behaviour of 
solids. As a result, calculations relying on (TIo^ render implausibly long times of tidal desp inning 
- see, for example the discussion and references in Castillo- Rogez et al. (2011). This tells us 
that in realistic settings the MacDonald-style approach f TH|) based on ( T15]) is inferior, compared 
to the Darwin and Kaula method ([2]), which may, in principle, be applied to any rheology. 
Despite this, the MacDonald torque remains a convenient toy model, capable of furnishing 
results which are qualitatively acceptable over not too long timescales (Hut 1981, Dobrovolskis 
2007). 



5 The MacDonald torque and the ensuing dynamical 
model by Goldreich. The case of an oblate body 

In this section, we shall trace how the MacDonald theory of bodily tides yields the model 
of near- resonance spin dynamics by Goldreich (1966). Then we shall recall an oversight in 
the MacDonald theory, and shall demonstrate that correction of that oversight brings a minor 
alteration into Goldreich's model of spin evolution. 

The goal of this section is limited to consideration of the tidal torque solely. So we assume 
that the body is oblate, and the triaxiality-caused torque does not show up. Appropriate 
for spinning gaseous objects, this treatment indicates that the tidal torque stays finite for 
synchronous rotation and vanishes at an angular velocity slightly faster than synchronous. 

In Section E] below we shall address the more general setting appropriate to telluric bodies, 
with triaxiality included. While some authors (e.g., Heller et al. 2011) ignore the triaxiality- 
caused torque in their treatment of solid planets, it turns out that inclusion of the permanent- 
figure torque renders important physical consequences and changes the picture completely, as 
will be seen in Section |6j 



5.1 The MacDonald torque 

We shall restrict ourselves to the case of the tidally perturbed secondary coinciding with the 
tide-raising one, so M sec = M* ec , and all the orbital variables are identical to their counterparts 
with an asterisk. 

Differentiating ( |14p with respect to A , and then setting A = A* , we obtain the following 
expression for the polar component of the torque, for low i : 

5 

T z = \GM 2 sec k 2 — - 3 sin (2 (z> - 6) At) +0{i 2 /Q) 

2 r (t) r* (t-At) v J 



(21) 

G M 2 sec k 2 \ sin ( 2 (z> - 6) At ) + 0(i 2 /Q) + 0(enAt/Q) , 



5 

3 „ , r o , R 



where the error 0(enAt/Q) emerges when we identify the lagging distance r*(t — At) with 
r*(t) = r(t) , as explained in footnote [5] in the preceding section. Referring the Reader to 
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Efroimsky & Williams (2009) for this and other technicalities, we would mention that (12 ip 
is equivalent to the Darwin torque only under the condition that the rheological model (Tl5l - 
I20aj) is accepted. For potentials, employment of model ( JT5l - I20ap enables one to wrap up the 
infinite series ([2]) into the elegant finite form For torques (truncated to / = 2 only), similar 
wrapping of the appropriate series is available within the said model. 

In the preceding section, we explained the geometric meaning of the longitudinal lag fTl7|) : 
its absolute value is the double of the geometric angle separating the directions to the bulge 
and the secondary as seen from the primary's centre. 

If we define a quantity "Q" via f l20ap . the MacDonald torque will look: 



5 

% = § G Ml c k 2 — - 3 sine + 0(z 2 /Q) 

1 r (t) r* (t-At) 



5 

| G M s 2 ec sgn(z> -0) + 0(i 2 /Q) + 0(enAt/Q) . 



(22) 



For a nonzero eccentricity, the quantity "Q" should not be interpreted as an instantaneous 
quality factor, because it is not guaranteed to interconnect the peak work or peak energy and 
the one-cycle energy loss in a manner appropriate to a quality factor - see footnote [8] in the 
previous section. Therefore a more reasonable and practical way of writing the MacDonald 
torque (12TI) would be through using (120b|) or (ITTj) : 



5 

% = ~ GM 2 sec k 2 ^ At2{v-9) + 0(i 2 /Q) + 0{enAt/Q) . (23) 

As we saw above, the MacDonald model is self-consistent only for At being a frequency- 
independent constant. However the factor v — 9 showing up in (123]) varies over a cycle, for 
which reason the torque needs averaging. This averaging is carried out in Appendix [Bj In the 
vicinity of the 1 : 1 resonance (for 9 close to n ) , the sign of v — 9 changes twice over a 
cycle, which makes the averaging procedure nontrivial. This situation is to be addressed in 
subsections 15.21 and 15.31 



5.2 Goldreich (1966): treatment based on the MacDonald torque 

In this subsection, we shall briefly recall a presently conventional method pioneered almost half 
a century ago by Goldreich (1966). Based on the MacDonald tide theory, this method has 
inherited both its simplicity and its flaws. 

The 1:1 resonance takes place when the spin rate 9 of the primary (the satellite) is equal 
to the mean motion n wherewith the secondary body (the planet) is apparently orbiting the 
primary. The formula (12"2"P for the MacDonald torque contains not the difference 9 — n but 
the difference 9 — v , for which reason the expression under the integral may twice change its 
sign in the course of one revolution. 

With aid of the formula 

5 

v = M + 2e sin.M + - e 2 sin2.M + 0(e 3 ) (24) 
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and under the assumption that e C n, we shape the difference of our concern into the form of 
9-u = (9-nj-2ne (cosM + ^ecos2Al^ + 0(e 3 ) 



cosAl + - e cos2Af 



or, equivalent ly, 



6 - v 



5 71 

cosAl + - e cos2Al h 0(e 2 

4 2ne 



(26) 



2ne 
Here 

rj = 9- M- u- n (27) 

is a slowly changing quantity, whose time-derivative r] = 6 — n becomes nil when the system 
goes through the resonance] 10 ! To impart the words "slowly changing" with a definite meaning, 
we assert that fj/n is of order e 2 - a claim to be justified a posteriori. 
Expression ( 1251) changes its sign at the points where 

cosM = - - e cos2Al + + 0(e 2 ) . (28) 

4 2 n e 

As all the terms on the right-hand side of (128!) are of order e at most, so must be the term 
on the left-hand side. Hence condition (|28|) is obeyed in the two points whose mean anomaly 
(and therefore also true anomaly) is close to ±ir/2: 

" = ± (| - ^) , (29) 
5 being of order e . From ( I29p and (I24p we obtain: 

sin<5 = cos^ = cos(.M + 2e sinA^ + 0(e 2 ) ) = cosA^ - 2e sin 2 Al + 0(e 2 ) , (30) 
which, in combination with (1281). entails: 



sin5 = --e cos27W + 2e sin 2 Al + 0(e 2 ) = e f - - - sin 2 Al ] + 0(e 2 ) 

4 2ne 2ne \4 2 J 

Insertion of (j29|) into ( 124|) also yields ± sinA^ = cos 5 + 0(e) . Recalling that 5 is of order 
e, we obtain: sin 2 Al = 1 — sin 2 5 + 0(e) = 1 + 0(e) . This enables us to rewrite (13~T1) as 

7 e + °( e2 ) • ( 32 ) 



2ne 4 



10 Goldreich (1966) defined this quantity simply as rj = 9 — M. , because he reckoned from from a 
fixed perihelion direction. Wc however reference our 9 from a direction fixed in space. (It is the same direction 
wherefrom the node is reckoned.) This convention originates from our definition of 9 as the sidereal angle - 
it is in this capacity that 9 was introduced back in equations ((3] - HJ . As we are not considering the nodal or 
apsidal precession, our subsequent formulae containing 9 will be equivalent to those ensuing from Goldreich's 
definition of 9 . 
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Before finding the rate fj = 9 — n at which the resonance is traversed, let us enquire 
if perhaps it could be simply put nil, the satellite being permanently kept in the resonance. 
The answer is negative, because for a vanishing 9 — n the difference 9 — v emerging in (125]) 
becomes a varying quantity of an alternating sign, and so becomes the torque. On general 
grounds, one should not expect that the average torque becomes nil for a vanishing rj , though 
it may vanish for some finite value of f] . 

To find this value of f] , many authors (Goldreich 1966, eqn. 15; Kaula 1968, eqn 4.5.29; 
Murray & Dermott 1999, eqn. 5.11) simply integrated the MacDonald torque (122]) . assuming 
the quality factor Q constant, and thus keeping it outside the integral: 



( Cioldreich) 



3GM* ec k 2 R 


1 




Ana 2 Q 


(1 - e 2 


) 1/2 J 


3GM* ec k 2 R 


2 




Ait a 2 Q 


(1 - e 2 


)l/ 2 


3 G M 2 sec k 2 R 5 


ptt/2-6 




2 7T a 6 Q [j 


'a ' J 


tt/2-5 



-2n 



R 



sgn( 



i>) dv 



r 



— t dv 

4 



— dv 



1 +4e cosv) dv + 0(e 2 



(33) 



3 G M spr k 2 R 5 



(4e cos 5 - 5) + 0(e 2 ) 



7r a 6 Q 

In a trapping situation, the average torque vanishes. So (134]) entails: 

<5 = 4ecos5 = 4e + 0(e 2 ) . 

By combining the latter with (132]) . the afore quoted authors arrived at 
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stall 



n e 



(34) 



(35) 



(36) 



an expression repeated later in papers and textbooks (e.g., eqn. 4.5.37 in Kaula 1968, or eqn. 
5.14 in Murray & Dermott 1999). This result however needs to be corrected, because the 
MacDonald approach is incompatible with the frequency- independence of Q assumed in (1341) . 

Be mindful that ( ]36]) . as well as its corrected version (144]) to be derived below, indicate that 
rj is of order e 2 . This justifies our assertion made in the paragraph after formula (|2T]) . 



5.3 The corrected MacDonald model 

To impart the MacDonald treatment consistently, one has to calculate the averaged torque, 
with the frequency- dependence of Q taken into consideration. As we saw in subsection H] the 
MacDonald approach fixes this dependence in a manner that can, with some reservations, be 
approximated with 



Q 



l 



(37) 
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or, in more general notations, 



Q = £ a x a , with a 



1 



(38) 



where the double instantaneous synodic frequency is given by (pi 

The form (138]) of the rheological law is more convenient, as it leaves us an opportunity to 
consider values of a different from —1 . For any value of a , the constant £ is an integral rhe- 
ological parameter, which has the dimension of time, and whose physical meaning is discussed 
in Efroimsky & Lainey (2007). It can be demonstrated that in the special case of a = — 1 the 
parameter £ coincides with the time lag At . For actual terrestrial bodies, a is different from 
— 1 , and the integral rheological parameter £ is related to the time lag in a more complicated 
manner {Ibid.). 

As demonstrated in Appendix [Bj insertion of (I38p and (1181) into (T2"2"j) . for a = —1 , or 
equivalently, direct employment of (1231) . entails the following expression for the orbit-averaged 
torque acting on a secondary: 



r z > 



- z 



9 A(e) - n Af(e) + 0(i 2 /Q) + 0(Q~ 3 ) + 0(enAt/Q) 



(39) 



3e 2 + V (l-e 2 ) 



where 
A(e) 
and 
Me) 

2 8 16 

while the factor Z is given by 



-9/2 



15 s 
1 + — 



105 



15 2 45 4 5 



(1-e 2 ) 



Z 



3 G M sec k 2 £ R 
R ^ 



3 n 2 Ml 



k 2 At R 



(M 3 



prim 



+ M SI 



27 2 

1 + Y e 



SnM sec k 2 



0(e 6 



573 



(40) 



e 4 + 0(e e 



(41) 



5 

R n 



Q (M prim + M sec ) a 3 x 



(42) 



M, 



prim 



and M sec being the masses of the primary and the secondary. 



Be mindful that the right-hand side of 



contains a multiplier ^ 



n 



- which is 
2| 6 - v | 

missing in the despinning formula employed by Correia & Laskar (2004, 2009). This happened 
because in Ibid, the quality factor was introduced as l/(nAt) and not as l/(x At) - see the 
line after formula (9) in Correia & Laskar (2009). In reality, the quality factor Q must, of 
course, be a function of the forcing frequency x (which happens to coincide with the mean 
motion n in the 3:2 and 1:2 resonances but differs from n outside these). 

The quality factor Q being (within this model) inversely proportional to x i the presence 

of the ^ factor in fT42|) makes the overall factor Z a frequency-independent constant. 

X 



Rewriting ( 1391) as 



%) 



Z 



Z 



15 9 105 

9 1 + — e 2 + 

2 4 

15 2 105 
V 1 + TT e + —T 



27 2 
n H 



573 



0(e 6 ) + 0(i 2 /Q) + Q(enAt/Q) 



6ne* 



121 . 

1 + e 2 
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+ 0(e 6 ) + 0(i 2 /Q) + 0(enAt/Q) , (43) 



11 Recall that, for 



-1 , the rheological parameter £ is simply the time lag At. 
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we see that it vanishes when the rate of change of r] = 6 — M. — uj — Vt accepts the value 
tun = 6ne 2 + | n e 4 + 0(e 6 ) + 0(i 2 /Q) + 0(enAt/Q) . (44) 

o 

For f] larger or smaller than rj staU , the average torque (1431) is nonzero and impels rj to evolve 
towards the stall value (j4*4"j) . 

On the right-hand side of (jSj), the leading-order term contains a numerical factor of 6, 
as different from the factor 19/2 showing up in fl36|) . The necessity to change 19/2 to 6 in 
the e 2 term was pointed out by Rodriguez, Ferraz-Mello & Hussmann (2008, eqn. 2.4), who 
had arrived at this conclusion through some different considerations (which, too, were based 
on the frequency- dependence (1551) ). The same result can be obtained within the Darwin-Kaula 
approach, provided the scaling law (1551) is employed (Efroimsky 2O12a,b).0 

Although correction of the oversight in the MacDonald torque renders a number different 
from the one furnished by Goldreich's development (6 instead of 19/2), qualitatively the princi- 
pal conclusion by Goldreich (1966) remains unchanged: when an oblate body's spin is evolving 
toward the resonance, vanishing of the average tidal torque entails spin slightly faster than 
resonant, a so-called pseudo synchronous rotation. It however should be strongly emphasised 
that the possibility of pseudosynchronous rotation hinges upon the dissipation model employed. 
Makarov and Efroimsky (2012) have found that a more realistic tidal dissipation model than 
the corrected MacDonald torque makes pseudosynchronous rotation impossible. 

6 Evolution of rotation near the 1:1 resonance. 
The case of a triaxial body 

In distinction from a gaseous or liquid body, a solid body would be expected to have a permanent 
figure in addition to the tidal distortion discussed so far. Goldreich (1966) demonstrated that 
this permanent figure plays an important role in determining if a primary, whose rotation is 
being slowed down or sped up by the tidal torque caused by a secondary, can be captured into 
the synchronous rotation state. This section follows Goldreich's derivation, but substitutes 
the MacDonald tidal torque with its corrected version, and also uses a more general mass 
expression. 

6.1 Rotating primary subject to a triaxiality-caused torque. 
The equation of motion and the first integral. 

Consider a triaxial primary body, which has its principal moments of inertia ordered as A < 
B < C , and which is rotating about the maximal-inertia axis associated with moment C . 
About the primary, a secondary body describes a near-equatorial orbit (so its inclination on the 
primary's equator, i, may be neglected). The secondary exerts on the primary two torques. 
One being tidal, the other is triaxiality-caused, i.e., generated by the existence of the permanent 
figure of the primary. Its component acting on the primary about the maximal- inertia axis is 

T tnax = \{B- A)°^ sin2A , (45) 

12 Implicitly, this result is present also in Correia et al. (2011, eqn 20), Laskar & Correia (2003, eqn 9), and 
in Hut (1981). The earliest implicit occurrence of this result was in Goldreich & Peale (1966, eqn 24). 



15 



the longitude A being furnished by formula (flOj) . In that formula, the sidereal angle 9 is 
reckoned from a reference direction in space to the principal axis associated with moment A . 
The acceleration of the sidereal angle then obeys 

C6--{B - A) sin2A = . (46) 

In neglect of the nodal and apsidal precession, as well as of , definition fl27|) yields: 

fj = 9 - n . (47) 

Ignoring changes in the mean motion0 we write: 

= V , (48) 

so the first term in (]46j) becomes simply Crj . 

To process the second term in (]46p . we would compare ffTUl) with (|27|) : 

A = - e + n + u + v + 0{i 2 ) = (-6 + n + uj + M) + {v - M) + 0(« 2 ) 



= - 77 + {y - M) + 0(z 2 ) (49) 

In neglect of the inclination, the following approximation is acceptable in the vicinity of the 
1:1 resonance 

(r- 3 sin2A) = - G 2 oo(e) a -3 sin 2r? , (50) 
where (...) signifies orbital averaging, while the eccentricity function can be approximated with 

G 2 oo(e) = 1 - \ e 2 + 0(e 4 ) . (51) 

This way, omitting 0(i 2 ) in (T49l and substituting sin2A with its average in fj46l) . we transform 
the second term in (T4"6"]) to: 

+ l(B-A) G 200 (e) sin 277 • (52) 

2 a 6 

13 The caveat about three neglected items implies that our n is the osculating mean motion n(t) = y M/ a (0 3 i 
and that we extend this definition to perturbed settings. The so-defined mean anomaly evolves in time as 
M = Mo{t) + f n(t) dt , whence M = Mo + n{t) . 

The said caveat becomes redundant when n is defined as the apparent mean motion, i.e., either as the mean- 
anomaly rate dAi/dt or as the mean-longitude rate dL/dt = dfl/dt + dui/dt + dA4/dt (Williams et al. 2001). 
While the first-order perturbations of a(t) and of the osculating mean motion y/ [i/a(t) 3 include no secular 
terms, such terms are often contained in the epoch terms Cl , w , and A4q . This produces the difference between 
the apparent mean motion defined as dL/dt (or as dAi/dt) and the osculating mean motion yj fi/a(t) 3 . 

14 Our neglect of evolution of the mean motion is acceptable, because orbital acceleration h is normally much 
smaller than the spin acceleration 9 . Indeed, for torques arising from the primary, n/9 is of the same order as 
the ratio of C to the orbital moment of inertia. 

15 Pioneered by Goldreich & Peale (1966), the approximation is explained in more detail by Murray & Dermott 
(1999) whose treatment omits terms of the order of e 3 and higher (see formulae 5.59 - 5.60 in Ibid.). To justify 
the omission, recall that before averaging the expansion of r -3 sin 2 A includes a series of terms with different 
arguments of the form sin(2?7 + qM) , where q is an integer. Since we are interested in dynamics in the vicinity 
of the 1:1 resonance, where r\ is a slow variable, then only the sin(2ry) term remains after the average. The 
averaged-out terms are of the order of e and higher powers including the e 3 terms. 
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While Goldreich (1966) assumed that the mass of the secondary is much larger than the mass 
of the primary, we do not impose this restriction. Combining Kepler's third law, G(M sec + 
M prim )/a 3 = n 2 , with formulae ( )46|) . (|48|) . and ( )52|) . we finally arrive at 

an equation describing the evolution of the rotation angle 77 . As pointed out by Goldreich 
(1966), this equation is equivalent to the one describing a simple pendulum. Indeed, in terms 
of 8 = 2?i , equation ( 1531) becomes 8 + y?., sin 8 = , with a constant positive vf ., 

'^■lib—max r — 1 '^1%0—max 
lCjl 



and with 8 playing the role of the pendulum angle. 

It follows directly from (153]) that the spin acceleration 7) vanishes if 77 assumes the values of 
or 71 . As can be seen from the pendulum analogy, the initial conditions (77 , rj ) t=t ^ = (0,0) 



as well as the conditions (77,?)) t=tQ = (it , 0) , correspond to the situation where the pendulum 

comes to a stall in the lower point (so 8 = when 8 = or 2tt). Under such initial conditions, 
77 stays or tt all the time, which implies a uniform synchronous rotation of the secondary 
about the primary. 

Multiplication of equation f )53|) by 7) , with subsequent integration over time t , gives the 
first integral of motion, 

\^^-\(B-A) M J'- Mwm * GMe) cos 2r/ — E , (54) 

whose value depends on the initial conditions. 

The period of the variable rj is given by a quadrupled integral over a quarter-cycle of rj : 

r r l~ r lmax fin 

P = 4 / f . (55) 

For circulation, use 77 max = tt/2 . In the case of libration about rj = , the value of 77 ma;c is less 
than 7r/2 and can be expressed via E by setting f] = in equation (I54p . To get an explicit 
expression of P , one should first express rj via E using f )5^|) and taking the positive root, and 
then should plug the so-obtained expression for fj into flo3]) . Also be mindful that libration 
about tt can be converted to the same integral by adding or subtracting tt from the variable 
of integration. 

In what follows, averaging over the period P will be denoted by (...) P , the subscript 
serving to distinguish the operation from averaging over the orbit employed in Section [5] and 
in formula ( |50|) . 



6.2 Parallels with pendulum. The auxiliary quantity W 

Goldreich (1966) noted the similarity of the differential equations f)53l- [5^|) to the classical pen- 
dulum problem. If the initial conditions place the body outside of the 1:1 spin-orbit resonance, 
then 77 circulates with a forced oscillation in rotation that depends on the mean value (7)) p 

16 In subsection l6.5l bclow. we write down the expression for the libration frequency Xub-max anc ^ a ^ so explain 
the reason why we equip it with such a subscript - see formulae ([55)1 and 

17 According to ()53l) . the spin acceleration 7) vanishes also for r\ = ±7r/2 , which is the upper point of the 
pendulum. 
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of the 77 frequency. If however the initial conditions place the body's spin within a sufficiently 
close proximity of the 1:1 resonance (the "resonance region"), then 77 will librate about or 
7r . This will be a free physical libration with an amplitude smaller than n/2 . Inside the reso- 
nance region, the initial conditions establish the free-libration amplitude and phase. However 
the restricted nature of the motion will keep the mean value of 77 constant: it will be either 
or ?T. For the same reason, (r]) p will remain zero in the libration regime. In contrast to 
this, outside of the resonance region the mean rate of circulation (r)) p will possess a value 
determined by the initial conditions on rj and fj . 

Employing f l5"4"|) to express rj via E , and plugging this expression into ( 15"5|) . Goldreich (1966) 
demonstrated that the period P can be expressed via E , inside or outside the resonance region, 
by a complete elliptic integral of the first kind. This is natural, as the expression of 77 through 
E involves a square root. 

Just as in the pendulum case, there exists a critical E for which the period diverges. This 
is the boundary E = Eb separating circulation from libration. To locate the boundary, one 
should set simultaneously 77 = 7r/2 and 77 = in formula ( I54p : 

E b = 3 -(B-A) - Msec G 200 (e) n 2 . (56) 

^ -'"sec *T J-'lprim 

Mathematically, the emergence of a logarithmic singularity in P at E — Eb can be observed 
from formulae (9 - 10) in Ibid. Physically, this situation resembles the slowing-down of a 
pendulum at the circulation/libration border. In our problem, though, this division corresponds 

tO 7] = ±7r/2. 

By setting simultaneously 77 = and 77 = in ( 154]) . we obtain the minimal value that the 
integral E can assume: 

minE = - E b . (57) 

The values E > E b correspond to circulation in either direction, with oscillations of the 
rotation rate. The values falling within the interval E b > E > — E b give libration. The 
minimum value E = — E b corresponds to synchronous rotation without libration. 

The quantity 77 being slowly varying (compared to the mean motion n), the period P , 
with which 77 is changing, naturally turns out to be much longer than the orbital period, both 
outside and inside the 1:1 resonance: 

P » — . (58) 
n 

This can be appreciated from the evident equalities 

„ 2 7T 

outside of the 1:1 resonance , (59a) 



2 71 

P = inside the 1:1 resonance . (59b) 

In (I59bp . Xub denotes the libration frequency which too is much smaller than the mean motion 
n , as we shall see shortly. 
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A useful quantity defined by Goldreich (1966) was 



W = P (rf) p = / rf dt = 4 / fjdrj . (60) 

J0 Jri=0 

Similarly to (IB3|) . r\ max = rc/2 for circulation, and r] max < n/2 for libration. Just as P , so 
W can be expressed through E by complete elliptic integrals. For circulating (rj) p , Goldreich 
described P and W as applying to one oscillation, but they describe one circulation of rj with 
two oscillations. For librations, P and W describe one libration cycle. 

One more quantity of use in this problem will be the mean square variation of fj about 
{rj) p , given by (if) p — (rj) 2 p . Be mindful that the notation (...) p is employed to indicate 
averaging over a circulation cycle as well as that over a libration cycle. 

When the values of the mean motion n , the mass factor M sec / (M sec + M prim ) , and the 
ratio (B — A)/C are given, and the initial conditions on 77 and r) are set, equation (1541) 
furnishes the value of E/C . The comparison of this value with Eb/C distinguishes circulation 
from libration. With the values of (B — A)/C and E/C known, the complete elliptic integrals 
can be evaluated and the values of the quantities P , W , (rj) p , and (rf) p can be found. It 
should also be possible to reverse this procedure. For example, the knowledge of the circulating 
value of (rj) p , along with the knowledge of n , the mass factor M sec j{M sec + M prim ) , and the 
ratio (B — A)/C , should allow E/C , W , and (rf) p to be computed. 



6.3 The tidal torque and the libration bias 

To account for tidal dissipation, Goldreich (1966) added the averaged (over an orbital period) 
tidal torque (134|) to the right-hand side of equation (1531) . We however shall employ the corrected 
average torque (|43|) instead of (|34|) . This will result in 



z 



L + y e 2 j - 6ne 2 M + ^ e 2 



0(e 6 ) + 0(i 2 /Q) + 0(enAt/Q) , (61) 



where the coefficient Z depends via (f42l) on the orbital variables, the tidal parameters, and 
the masses. 

As we already mentioned at the beginning of Section [5j ignoring the permanent-figure term 
would lead one to the conclusions that the tidal torque is finite for synchronous rotation and 
that it vanishes for a spin rate slightly higher than synchronous. However inclusion of the 
permanent-figure term into the picture alters the results radically. The constant term on the 
right-hand side of equation (IBTj) causes the mean value of 77 to be slightly larger than or tx . 
For small librations, this bias is 

= J^Wn A -H^ (> + f ^) G^) + ° <e6) + ° { ' VQ) + ° {mAt/Q> <62a) 



2 Z e 
(B-A)n 



M. 



sec 



1+^e 2 



(62b) 
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with Z calculated through (142|) . The value of the time lag At entering the expression (1421) 
for Z should be set appropriate to the orbital frequency (mean motion). 

As evident from (1621) . the bias comes into being due to the eccentric shape of the orbit. 
Bodies with permanent figures can achieve synchronous rotation because the bias in 77 gives 
birth to a permanent-figure torque that balances the constant dissipative torque. 

At this point, an important caveat will be in order. As we explained in Section HI the 
MacDonald model of tides becomes self-consistent only when the time delay At emerging in 
(1421) is set frequency- independent. This circumstance limits the precision of the MacDonald 
torque and of the Goldreich dynamical model based thereon, whenever the model gets employed 
to determine the timescales of spin evolution (Efroimsky & Lainey 2007). Nevertheless, for very 
slow evolution the model may be employed for obtaining qualitative estimates. In this case, At 
should be treated as a parameter that itself depends upon the forcing frequency in the material. 
For rotation outside the 1 : 1 spin-orbit lock, it would be a tolerable approximation to use 
At appropriate to the principal tidal frequency X2200 or ^° ^ ne double instantaneous synodic 
frequency (1181) . However inside the 1 : 1 resonance, At would correspond to the libration 
frequency \ lih which may be very different from the usual tidal frequencies for nonsynchronous 
rotation. This circumstance may change the value of At and therefore of Z noticeably. 

6.4 Tidal dissipation and the point of stall 

Be mindful that on the right-hand side of equation flBTj) we have the tidal torque averaged 
over the orbit (see Appendix B for details). Similarly, the second term on the left-hand side of 
(IBTT) is the permanent-figure torque averaged over the orbit - see expression Therefore 
equation (l6Tj) renders us the behaviour of rj over times longer than the orbital period. This 
is acceptable, because the orbital period is much shorter than the timescales of our interest. 
Specifically, it is much shorter than P , see equation (J5HD • 

Without dissipation, E is conserved during oscillations of 77 , while in the presence of weak 
dissipation E changes slowly. Indeed, multiplying both sides of (l6Tj) by 77 and making use of 
(1541) . we arrive at 

+ 0(e 6 ) + 0(i 2 /Q) + 0(enAt/Q) 

Thus we see that a dissipative tidal torque influences the spin, while causing changes also in 
the first integral E (which is not identical to the actual kinetic energy). Considering equation 
( j54|) and noting that the cosine term is periodic, we see that weak tidal dissipation will cause 
a change in rj 2 over time scales long compared to P . 

For nonsynchronous rotation, weak dissipation causes both a slow secular change and a 
small oscillation. The latter arises from the oscillating part of r) as 77 circulates or librates. 

Goldreich (1966) also averages over the period P . For libration, the so- averaged rate of 
change of E is 

(dE/dt) p = - (\ + y e 2 ^j + 0(e 6 ) + 0(i 2 /Q) + 0(enAt/Q) , (64) 

where the positive definite quantity W is introduced via f )60|) . The negative (dE/dt) p for 
libration means that the maximum fj 2 at cos2/7 = 1 in equation (l54"j) is decreasing and the 

18 For the first term, Crj , the caveat about orbital averaging is unimportant, because rj bears no dependence 
upon the mean anomaly. 
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free libration damps with time. Libration evolves toward synchronous rotation with rj = and 
rj equal to or to ir . 

The rj term causes the damping of the libration amplitude. When the amplitude is suffi- 
ciently small, its decrease obeys the exponential law exp(— D L t) with 

D ^i( 1 + f e2 ) • < 65 > 

As we emphasised in the paragraph after equation (ISTj) . inside the 1 : 1 resonance the time 
delay At showing up in the expression for Z should be appropriate to the libration frequency 
X Ub which may differ greatly from the usual tidal frequencies for nonsynchronous spin. This 
choice of At will influence the value of Z . 

For circulation, averaging of equation ( [63]) over one orbital period leads to 

(dE/dt) p = - | 



W 



15 , 

1 + Y e ' 
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+ 0(e 6 ) + 0[i 2 jQ) + 0{enAt/Q) (66) 

where we have used definition f)59ap for P and definition (160]) for W . 

While W = P(r] 2 ) p is positive definite, the second term on the right-hand side of equation 
(1661) can, for circulation, have either sign. In the case when its sign is positive, the expression 
( |66l) for ( dE jdt) p will vanish for W equal to 

W stall = 12 nn e 2 (l + ^) + 0(e 6 ) + 0(z 2 /Q) + 0(enAt/Q) . (67) 



No matter whether W begins its evolution with an initial value larger or smaller than W, 



stall ' 

it will never cross W stall . 

Another important value of W is the one corresponding to the boundary between libration 
and circulation, Wb . Below we shall obtain its value and shall explain that for Wb < W stall 
there exists a positive value of (77) p at which the evolution of a circulating (77) p should stall. 

6.5 The free-libration frequency 

Exploring librations, we start out with the small-amplitude case. Replacing sin 2rj with 2r\ in 
equation fl53l) . we write down the frequency for small librations of 77 : 



2tt 

^lib—max 



P„ 



B-A Msec 

3 ~C~ M sec + M~ G20o(6) 
B — A M st 



1/2 

n . (68) 



In many realistic situations, the condition ^ ri^ -jj se \/r 1 is fulfilled0 which 

ensures that Xm-max ^ n • Larger amplitudes increase the period P and render a smaller 
frequency x lih , so expression flS"8l gives the smallest librating period and the largest frequency 

Xlib-max = ma *Xub ■ (69) 



19 This condition is satisfied safely if the primary is a planet (like Mercury) or a large satellite (like our Moon). 
Its fulfilment is not guaranteed, though, for small satellites (like Phobos or Hyperion). 
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The linear 7) term in (|6T|) will alter the frequency, but for slow tide-caused evolution the 
correction will be small, so ( 168]) still will serve well as an approximation for the maximal 
frequency of libration. On all these grounds, we now accept that for both small-amplitude and 
large-amplitude librations the inequality 

X lib « n (70) 



holds. This justifies, a posteriori, the assertion made after (159b|) . 



In the literature on the physical libration of the Moon, the expression for the free-libration 
frequency is ubiquitous, though often without the mass factor. Versions of the expressions for 
Vbias an d the free-libration damping rate appeared in Williams et al. (2001). That paper, 
though, did not rely on the corrected MacDonald model (constant time delay), but simply used 
separate Qs for libration and orbital frequencies. 

6.6 The boundary between circulation and libration 

The boundary between circulation and the resonance zone corresponds to a value W = Wb . 
To find it, we combine (I54p with (156 p and arrive at 

I C1) 2 = E b (1 + cos2r/) , (71a) 



2 



which is the same as 



7)2 = 4 COS 2 T] . (71b) 



Insertion of the resulting expression for f] into (1601) entails: 



^ r /2 . _ [e, 



Wb = 8 \ j — I cos rj drj = 8 y — = An 



'0 




B - A M, 



sec 



C 2 oo(e) 



1/2 

, (72) 



comparison whereof to ( 1681) yields another expression for the boundary value: 

**b — 4 Xlib-max 

While W is continuous across the boundary, P has a logarithmic singularity, as was mentioned 
in the paragraph after equation ( 1561) . As P diverges, the evolution rate of (dE/dt) p , given 
by ( 1661) , vanishes at the boundary. Nonetheless a small perturbation allows the boundary to 
be crossed - for more on this, see the two paragraphs after formula (22) in Goldreich (1966). 

6.7 Three regimes 

If W stall < Wb , then there is no stall point in the region of circulation (W& < W). So W 
evolves towards Wb , while (77) p of either sign evolves towards zero at the libration/circulation 
boundary. Figure la illustrates the evolution of (r)) p . Goldreich (1966) comments that the 
boundary will be crossed, the free libration will damp, and the rotation will evolve toward the 



22 



synchronous state. The synchronous state has a zero n, with 77 biased slightly offp of either 
or 7r . The figure does not show the libration region] 21 ! 

If Wb < W staU , then a stall point exists in the circulating region ( Wb < W). All in all, Figure 
lb summarises the following three cases: 



A. For an initially negative (77) , circulation is taking place. We have Wb < W , but the stall 
point is never located on the negative (77) p side of zero. The rate ( dE/dt ) p is negative, as 
can be seen from fl63l) . We then have a slow decrease of the three quantities!^ E — >■ Eb , 
W — > Wb , and (f)) p —> 0. This decrease takes the system to the circulation/libration 
boundary. Goldreich (1966) commented that it is ambiguous whether (77) p would cross 
zero and proceed to increase or whether the boundary would be crossed passing into 
the libration region followed by damping of the libration and by evolution toward the 
synchronous state. 

Be mindful that, while the first integral E is decreasing in the circulation case, the ac- 
tual kinetic energy of rotation is increasing. Indeed, when the negative (f)) p is evolving 
towards zero, the spin rate 9 is growing, and so is the rotational energy C9 2 /2 . 



B. For an initially positive (77) p lying between the circulation/libration boundary and the 
stall point, i.e., for Wb < W < W stall , the rate (dE/dt) p is positive, and E increases. 
The quantities W , (f] 2 ) p and (fj) p will evolve to larger values until the evolution stalls 
as W approaches W stall from below. 

As the positive i) is increasing, so are the spin rate 9 and the rotational energy^ 

C. For an initially positive (?)) p beyond the stall point we have: Wb < W stall < W. Then 
(dE/dt) p is negative and E is slowly decreasing! 2 ^! Likewise, W , (?) 2 )p and (77) p 
evolve to lower values until the evolution stalls as W approaches W stall from above. 

The decrease of the positive rj renders a decrease in the rotation rate 9 and thus leads 
to damping of the kinetic energy of rotation. 



20 The value it shows up because of the factor 2 accompanying 77 when this variable enters sin 2?? in the 
equation for the torque and cos 2rj in the expression E . The tiny bias off rj — or off i) = 7T , given by (l62l) , 
will emerge due to dissipation. 

21 A figure with W rather than (ry) p would show the libration region, along with the circulating region. On 
such a figure, though, positive and negative values of (rj) p would overlap, so that the one-sided nature of the 
stall would not be evident. 

22 Be mindful that E remains a constant over a cycle of P , except for a tiny amount of dissipation during 
a cycle. 

23 In the first integral ([54]) . the kinetic-energy- like part is given by Cr] 2 /2 = (C/2)(#-n) 2 = (C/2)(9 2 - 2n8 + 
n 2 ) . When the spin accelerates, the additional kinetic energy is borrowed from the orbital motion. Because 
of the dissipation, the overall spin + orbit energy must nevertheless decrease. (The orbital energy includes a 
kinetic part and the —GM/r potential term.) 

24 Just as in item A, here E stays virtually unchanged over P . 
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Calculation of the frequency at the stall point (f} staU ) P would require two steps: first, 
reversing the elliptic integral computation, starting with W sta ii] and second, deriving E/C , 
Pstaii and (f) stall ) p . 

At the stall point, formulae ( I59al) and (IBTfl) acquire the form of 

2 7T 

P stall = -J^ T aIld W ^all = P stall )p > ( 73 ) 

V I stall I P 



correspondingly. In combination with (l67j) . this entails: 

\ 'stall' P c„„2 ' - e 



6ne 2 1 + - , (74) 

whence we once again note that the frequency {f) sta u ) P must be positive. The comparison of 
^ ' stall w hh in the inequalities mentioned in the above items A, B, C is then equivalent 
to comparing (77^ ) p /{r) stall ) p with (2/tt) x^_ max or to comparing^ 3 tt e 2 (1 + e 2 /r6) with 
q B — A M Earth p / \ 1 1 ' 2 

<^ Ml Earth + ^Moon J 

A comment on equation (!74|) would be in order. Since oscillations of 77 result from the 
existence of the permanent triaxiality, evolution of 77 becomes smooth in the A = B limit. 
Averaging becomes unnecessary, so (Vs ta u)p becomes simply 77 2 „ , while (i) staU ) P becomes 
fj stall ■ This way, in the oblate-body case considered back in Section El equation ( 1741) acquires 
the form of 77 sta „ = Qne 2 (1 + e 2 /16) , which agrees with (j44j). 



6.8 Application to the Moon 



When Goldreich (1966) applied his inequality expressions to the Moon, he found that the stall 
point would have interrupted evolution of rotation from faster spin to synchronous rotation. 
Here we have repeated his study, though with the corrected average torque (l4"3l) instead of (134]) . 

For B 7A = 2.278 x 1(T 4 (Williams & Boggs 2012) , e = 0.0549 , and M£ , arf A = 

^ J-Vl Earth + M Moon 

0.98785 , the libration period turns out to be 38 times the orbital period, while W ata ii is 
10% larger than Wb . This renders a value of W s t a ii/Wb much smaller than the one found 
by Goldreich (1966), and the difference is mainly due to our use of the corrected average 
torque (j4"3|) . Despite the so-different value of W sta ii/Wb , despinning of the Moon would still be 
interrupted by a stall. In principle, a tidal spin-up scenario remains an option too, though this 
option does not look probable. 

Goldreich (1966) noted that the lunar orbital eccentricity is changing. To make the stall 
point disappear, i.e., to ensure that W sta u < Wb- the eccentricity of the Moon would need to 
be less than 95 % of its present value. While we lack data on the ancient evolution of the lunar 
eccentricity, the modern eccentricity rate of about 2 x 10" 11 yr _1 has been reliably determined 
through analysis of the Lunar Laser Ranging data (Williams et al. 2001, Williams & Boggs 
2009). For the measured eccentricity rate, the eccentricity would be small enough to prevent a 
stall prior to 1.4 x 10 8 yr ago. The Moon has clearly been a satellite of the Earth for billions 
of years, and the capture into the synchronous spin state should have occurred very early in its 
history. 
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The latter may also be expressed as comparison of 37re' 



(i + ^ 2 ) 



with 
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M Ea 



Earth 
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For the Moon as it exists today, the free libration in longitude has a 2.9 yr period. The 
damping time is four orders of magnitude longer (Williams et al. 2001), so evolution of this 
libration is slow. Despite the damping time being short compared to the lunar age, the Moon 
has a small free libration amplitude of 1.3" (Rambaux and Williams 2011). There has been 
geologically recent stimulation, probably due to resonance crossing (Eckhardt 1993). 

The ambiguity of evolution of the negative (rj) p past the circulation/libration boundary 
deserves a comment. For Mercury, Makarov (2012) finds that a more realistic tidal dissipation 
model than the corrected MacDonald torque strongly changes computations of the evolution of 
planetary spin rate near the 3 : 2 and higher spin-orbit resonances. 

7 Conclusions 

In the article thus far, we have provided a detailed explanation of how the empirical Mac- 
Donald model can be derived from a more accurate and comprehensive Darwin-Kaula theory 
of bodily tides. We have demonstrated that the derivation hinges on a key assertion that 
the quality factor Q of the primary should be inversely proportional to the tidal frequency. 
This crucial circumstance was missed by MacDonald (1964), who made his theory inherently 
self-contradictory by setting the quality factor to be a frequency- independent constant. 

We have corrected this oversight in the MacDonald approach, and have developed an appro- 
priate correction to Goldreich's model of spin dynamics and evolution near the 1:1 spin-orbit 
resonance. Although we got different numbers, qualitatively the main conclusion by Goldreich 
(1966) stays unaltered: when an oblate body's spin is evolving toward the resonance, vanishing 
of the average tidal torque still implies a pseudosynchronous rotation (rotation slightly faster 
than resonant), while synchronicity requires a small compensating torque. For a triaxial body, 
the picture gets more complex due to the emergence of a triaxiality-caused torque. (While the 
oblate case is appropriate for spinning gaseous or liquid planets and moons, the triaxial case 
applies to rocky objects.) 

Goldreich (1966) linked the possible trapping of a body in the synchronous state during 
its tidal evolution of rotation to its triaxiality. In light of the correction of expression ( 13"6"|) to 
(I44j) . that limiting condition between the triaxiality and e 2 must change. After capture into 
the synchronous state, the tidal torque is compensated by a triaxiality torque by aligning the 
principal axis associated with the smallest moment of inertia slightly off of the mean direction 
to the external body. A constant is thereby introduced into the physical libration in longitude. 

Setting the tidal quality factor to scale as inverse frequency is incompatible with the actual 
dissipative properties of realistic mantles and crusts. Nevertheless, after the afore-explained 
correction is implemented, both the corrected MacDonald description of tides and the dynamical 
theory based thereon remain valuable toy models capable of providing a good qualitative handle 
on tidal dynamics over not too long timescales. Specifically, this approach renders a simple 
qualitative description of the interplay between the tidal torque and the triaxiality-caused 
torque exerted on a body near the 1:1 spin-orbit resonance. 

It should also be remembered that the stability of pseudosynchronous rotation hinges upon 
the dissipation model employed. Makarov and Efroimsky (2012) have found that a more real- 
istic tidal dissipation model than the corrected MacDonald torque makes pseudosynchronous 
rotation unstable. Finally, pseudosynchronism becomes impossible when the triaxiality is too 
large (see subsection 6.7, specifically footnote 24). 
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Appendices. 

A The tidal-torque vector 

and its components in spherical coordinates 

A secondary with spherical coordinates r* = (r*, A*, </>*) and mass M* ec raises a tidal bulge 
on the primary. The gravitational attraction between the tidal bulge and a secondary at 
r — (r , A , (p) with mass M sec causes equal but opposite torques on the primary and the 
secondary. For the external tidal potential U(r) , the torque components depend on the partial 
derivatives of the potential U along great circle arcs. 

To calculate these expressions, let us recall some basics. The torque T wherewith the 
primary is acting on the secondary is given by the cross-product 

f = fxf , (75) 
J 7 being the tidal force exerted by the primary on the secondary. This force is given by 

f M M ( — e - — !— — e - — e ^ (76) 

sec dr sec \ dr r r cos0 OX A r d<f> ^ J 

e r , e A , being the unit vectors of a spherical coordinate system associated with the primary's 
equator and corotating with it@ Insertion of f J76|) into f J75|) results in 

f = -M, ec ( e r+ ^e J+ ^f «,) , (77) 

The torque r wherewith the secondary is acting on the primary will be the negative of (1771) : 

? = M «« (°^ + % *» + sbf «*) < 78 > 



26 Were we using the polar angle (p = ir/2 — <j> instead of the latitude <f>, the right-handed triple 
of unit vectors would be: e r , , ("radial - south - east"), while the gradient would read as 

dU(f) dU a_i_1 3U 1 dU a 

— W~ _ ~W G r + r Sip e <P + f im^ ~$X e >- ■ 

As we are employing the latitude, the right-handed triple changes to e r , e A , ("radial - east - north"), 

and the gradient becomes = ^ e r + i ^ e A - I |^ e . 



26 



Its east component, the one aimed along e A , is parallel to the primary's equatorial plane. 
The north component, aimed along e,, is tangent to the meridian. In our paper however we 
employ the projection of the torque vector onto the primary's spin axis, i.e., the component t z 
orthogonal to the equator plane. This component is cos0 times the north component: 

dU 

% = M sec — . (79) 

This torque component, equation ([6]), slows down the rotation rate of the primary. The decel- 
erating torque acting on the secondary has an opposite sign and is rendered by (131) . 

The east component and the projection of the north component onto the equator plane act 
to change the orientation of the primary's spin axis, a subject we shall not pursue in this paper. 



B Calculation of the average torque within the corrected 
MacDonald model 

To calculate the orbital average of the tidal torque acting on a librating secondary obeying the 
corrected MacDonald tidal model, substitute ( J38l) and (fTBl into ( l22l) . and then choose a = — 1 
(and recall that, for a = — 1 , the integral rheological parameter £ is simply the time lag: 
£ = At). An equivalent option would be to employ ( 1231) directly. This will lead us to the 
following expression for the torque: 

5 

% = - ~ GM 2 ec k 2 %£ X sgn(fl + 0(z 2 /Q) + 0{enAt/Q) + 0(Q~ 3 ) 

5 

= - | GM 2 sec k 2 \ 2£ \9-v\ sgn(# - v) + 0(i 2 /Q) + 0{enAt/Q) + 0(Q~ 3 ) 

5 

= - 3 G M 2 ec k 2 ^ £{6-v) + 0(i 2 /Q) + 0(enAt/Q) + 0(Q" 3 ) , (80) 
r 

and for its average over one orbiting cycle: 

2 6 

(%)= - ~~~~^r hS + 0(i 2 /Q) + 0(Q~ 3 ) + 0(enAt/Q) = 

2 6 2 6 

- 3GMs : c 9 ( \ ) + 3GM Z C hS ( y% ) + O^/Q) + 0(Q~ 3 ) + 0{enAt/Q) (81a) 
K r° R r b 

3GM* ec k 2 £ ■ R & , 2 ,_ 9/2 1 [ 2 - 4 
= Q-{l-e) -j Q (1 + ecos,)^ 



+ 3 G M sec fc ^ n ^(!_ e ^ 6 i, 1* (l + e cos t/) 6 di/ + 0{i 2 /Q) + 0(g" 3 ) + 0(enAt/Q) . (81b) 
Evaluation of the integrals is straightforward and entails fl39l - H2l) . 
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Figure 1: Three possible scenarios of evolution of (rj) p . Arrows on the top diagramme illustrate the 
evolution of circulating (f]) p toward the libration boundary for W t „ < . Arrows on the left of 
the bottom diagramme depict the evolution of negative {rj) P toward the libration boundary. Arrows 
in the midst and on the right of the bottom diagramme show the evolution of positive {rj) p toward 
the stall point for Wb < W tall . 
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